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We show that non-abelian potentials acting on ultracold gases with two hyperfine levels 
can give rise to ground states with non-abelian excitations. We consider a realistic gauge 
potential for which the Landau levels can be exactly determined: the non-abelian part of the 
vector potential makes the Landau levels non-degenerate. In the presence of strong repulsive 
interactions, deformed LaughUn ground states occur in general. However, at the degeneracy 
points of the Landau levels, non-abelian quantum Hall states appear: these ground states, 
including deformed Moore-Read states (characterized by Ising anyons as quasi-holes), are 
studied for both fermionic and bosonic gases. 



The quest for states of matter having non-abelian excitations has been very intense in the past 
two decades, motivated by interest for the topological properties of such correlated states yj, |2[] 
and by their relevance for the implementation of topological quantum computation schemes ||20. 
Quantum Hall systems provide a major arena in which investigating non-abelian anyons, and a 
huge amount of work has been devoted to characterize which quantum Hall states are non-abelian 
and to propose (mainly interferometric) experiments to test the non-abelian nature of such states 
[|2|]. Although indirect evidences has been found for the v = 5/2 quantum Hall state [jSi |4|], a 
completely unambiguous signature of its non-abelian nature is still lacking and the experimental 
manipulation of anyons in such systems stands so far as an important open challenge M\- 

A significant and complementary issue consists in finding other experimental systems suitable 
to simulate the quantum Hall physics and realize non-abelian anyons. In this respect, ultracold 
atomic systems provide a natural candidate [|5i l6|], due to the possibility of using them as sim- 



ulators of many -body systems [|7|]. Two essential (but in general not sufficient) ingredients are 



available in ultracold atomic systems. First, one can simulate artificial magnetic fields by using 
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rotating traps iQ] or spatially dependent optical couplings between internal states of the atoms 
Correspondingly, the possibility to have Laughlin states (with abelian excitations) using strongly 
interacting, rapidly rotating two-dimensional ultracold Bose gases has been discussed in the liter- 
ature and it has been showed as well that incompressible vortex liquid states obtained for rotating 
ultracold bosons can be described by the Moore-Read wavefunction [5] (see the review [SI]). The 
second ingredient, useful in the perspective to have ground states with non-abelian excitations, 
is given by the possibility of trapping, both for ultracold bosons and fermions, two (or more) 
hyperfine levels of the same atomic species. 

By combining the two ingredients, one may think to realize a system with two hyperfine levels, 
coupled between them and each one feeling a different artificial effective magnetic field. Several 
possible schemes have already been proposed in literature to realize such non-abelian gauge po- 
tentials, acting on the two (or more) components of an ultracold atomic gas: in [jlO] laser fields 
couple three internal states with a fourth (the so-called tripod scheme), in such a way that for 
two degenerate dark states one can have a truly non-abelian gaug e potential. The extension of 
this scheme to the tetrapod setup has been recently discussed Another proposed scheme is 
based on laser assisted tunneling depending on the hyperfine levels to obtain U (2) vector poten- 
tials acting on ultracold atoms in optical lattices [|l2n . while effective non-abelian gauge potentials 
in cavity QED models has been recently addressed in [ isj]. Several properties of ultracold atoms 
in artificial non-abeli an gaug e potentials, including Landau levels and dynamical regimes, have 
been recently studied il4l-l2ll]. 

Motivated by these proposals, the natural question is whether one can find non-abelian anyons 
in such systems. Generally a non-abelian gauge potential in presence of strong interactions does 
not guarantee the non-abelianity of the excitations; rather, deformed Laughlin states (characterized 
by abelian quasi-holes) appear. Recently, different kinds of anomalous quantum Hall effect for 
ultracold atomic gases in artificial gauge potential have been addressed [Il9l - l2lll : in particular, in 
Um it has been shown that a ultracold Fermi gas in a lattice subjected to a non-abelian gauge field 
can present non-chiral and anomalous quantum Hall effects. 

In this paper we consider a non-abelian potential in a symmetric gauge acting on a two- 
component ultracold gas: this realistic potential (realizable with a suitable choice of the Rabi 
frequencies in a tripod scheme) has the advantage that Landau levels with high degeneracy are 
present and can be exactly determined. The two Landau levels become non-degenerate for a non- 
vanishing strength of the non-abelian potential: we show that, as expected, deformed Laughlin 
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states arise. However, at the points in which there is degeneracy of the Landau levels, ground 
states with non-abelian excitations do appear and explicit analytical results can be obtained. These 
ground states, including deformed Moore-Read states (presenting Ising anyons), are studied for 
both fermionic and bosonic gases. 

Single-particle Hamiltonian: We first consider the two-dimensional non-interacting Hamilto- 
nian for atoms in two different hyperfine levels (hereafter denoted by |t) and \\)) in a non-abelian 
gauge potential A: 

H={p, + A^f + {py + Ayf = Ha + Hna- (1) 

In the Hamiltonian ([T]) we explicitly introduced the (abelian) term proportional to the identity. Ha, 
and the (non-abelian) off-diagonal term Hna (the units are taken to have mass 1/2). 
We choose the following vector potential A: 

B B 

Ax = qa^ ~ ^ ^y = ^^y + y^' (2) 

where the a's are the 2x2 Pauli matrices. As we show in the Appendix, it is possible to engineer 
this vector potential in a system of rotating atoms having three quasi-degenerate ground states 
coupled with an excited state (the so-called tripod scheme): the vector potential A is in general 



determined by the choice of the parameters of the Rabi frequencies f2j fllOD. In the Appendix we 
show that the Rabi frequencies in a rotating trap can be arranged to give the vector potential 

In ^ the parameter q measures the strength of the non-abelian contribution and B is an artifi- 
cial, diagonal magnetic field orthogonal to the system. When q = 0, the usual Landau levels [|80 
are retrieved and are doubly degenerate. The total effective magnetic field is 

B = V X A + iAx A = B -2q'^a,. (3) 

One can see that the system is characterized by a translationally invariant Wilson loop jiol . 

To diagonalize the Hamiltonian we introduce the complex variables z = x — iy and z = 
X + iy, and we rewrite Ha and Hna as 

1 , / -id\ 

Ha = 2q^ + B + -dU, H,,a = q\ (4) 

4 yzd^ J 

where we introduced the creation and annihilation operators = Bz — Adz and d = Bz + Adz- 

The non-abelian part of the Hamiltonian Hna corresponds to a Jaynes-Cummings coupling be- 
tween subsequent Landau levels with different internal degree of freedom. Its spectrum is similar 
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to the one obtained in the relativistic case typical of the graphene systems and, in particular, Hna 
is analogous to the Hamiltonian obtained in ll20ll starting from the Dirac equation in an anisotropic 
regime. Therefore we can notice that it leaves a single family of uncoupled states, corresponding 
to the lowest Landau level, ipo \D, and, otherwise, its eigenstates take the form ^n-i It) i 




FIG. 1: Energies of the first 6 eigenstates as a function of q^/B. The n^^ Landau level is split by the non- 
abelian contribution of the external field in Xn+i (black solid) and xt (r^^ dashed); for g = one recovers 
the usual Landau levels, while for q^/B = 3 there is a degeneracy of the Landau levels (indicated by a 
circle). 



Using these wavefunctions we can diagonalize the whole single-particle Hamiltonian H ob- 
taining the following eigenvalues 



= 2Bn + 2q^ ± y^B^ + 8q^Bn 



(5) 



corresponding to the unnormalized eigenstates 



'X.r, 



(^B + 2qV2B^ T + 8q^Bn^ |t) + (^B - 2qV2B^ ± ^/B^~+8q^^ ipn \i) 



(6) 



where the ijjk are wavefunctions in the k^^ Landau level such that 2\^2Bn^pn = iSipn-i, due to the 
Jaynes-Cummings coupling. Therefore there is a correspondence between the usual eigenstates in 
a magnetic field and the states x'- every Landau level is split in two parts corresponding to the 
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states xt-i ^iid Xn ^^'^ for g — )■ their energy is ~ En-i ± 4g^n. The energy spectrum ^ is 
plotted in Fig HI 

The uncoupled state family corresponds to Xo energy, B + 2g^, is higher than the 

energy Ei of Xi^ which is the ground state family of the system for < 3B (the general case 
with > 3B will be analyzed in the following). To describe the states in the family Xi it is useful 
to introduce the operator Qi = c^^iCTx + ic/^ where forn = 1, 2 . . . we defined the constants: 



Qi allows to map uncoupled states in Xo into states in Xi so that we can describe every ground 
state in the form 



where P is a generic polynomial in z. 

Two-body interactions: We consider now a system of interacting atoms. We will treat both 
fermions and bosons: in order to obtain quantum Hall states, a strong repulsion between atoms 
of the same hyperfine level is needed. For fermions the intra-species repulsion is provided by the 
Pauli principle; for bosons, we assume an (intra-species) interaction between |t) atoms or between 
II) atoms of the form Hj = g J2i<j ^ i^i ~ ^i)' where (7 is a positive coupling constant such that 
g ^ B,q (the s-wave scattering lengths between |t) atoms and between ||) atoms are assumed 
equal). For the sake of simplicity, we assume as well that the inter-species interactions (i.e., 
between |t) and ||) atoms) is vanishing: in the following we will also comment on the stability of 
the obtained ground states with respect to a non-vanishing repulsive inter-species interaction. 

In order to find a ground state of the multiparticle Hamiltonian Ti = J^i' + must 
exploit the single-particle ground state degeneracy to obtain a wavefunction where all the particles 
lie in a superposition of states in Xi whose components |tt)ij ^nd ||i)jj vanish if Zi — )■ Zj. For 
fermions one finds (for < 3B) that, given a Laughlin wavefunction 



Ct,n 



B + 2qV2B^ + y^B^ + Sq^Bn 





(7) 



N 



A!r^ = n (-^ 



.,)™e-fS^N^ III ...I) 



(8) 



i<j 



with m odd, then the state 



N 



(m) 



(9) 



■N 



J 
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is a ground state of the Hamiltonian "H: every particle lies in a superposition of states Xi and the 
antisymmetric wavef unction causes the intra- species interaction energy to be zero (here and in the 
following Qn;j denotes the operator Qn applied to the particle j). For repulsive bosons, the form 
dH) holds with m even and greater than 2 (the states with m > 2 can be shown to be stable 
also under the presence of an inter-species repulsive interaction). 

The state \E'*^™^ can be seen as a deformation, due to the non-abelian potential, of the usual 
Laughlin states and it describes an incompressible fluid of spin-1/2 particles: its norm can be 
easily written in terms of the norm of the corresponding Laughlin state, showing that this state 
has, in the thermodynamical limit, a constant density. Also in the presence of quasi-holes in the 
Laughlin state 

N 

^tS = n ^1;^ - CO' (^^ - C2)' A!r^ (10) 



the norm is related to the corresponding quantum Hall state: the Berry phase due to the adiabatic 
exchange of the quasi-holes (^i and (^2 is the same of the one characterizing the corresponding 
quasi-holes in a simple Laughlin state. These excitations and their braiding statistics are therefore 
abelian [|22|] . 

Higher values of q^/B: The previous results can be extended to q'^/B > 3, where the 
ground state of the single-particle Hamiltonian is no longer Xi' one has that Xn>i '^he ground 
state for (2n — 1) < q^/B < (2n + 1) (see Figj2l). It is useful to introduce the operators 
Qn = c^,n d'^^'"~^^ax + c^,n d'^ ^ foT cvcry n so that the ground state wavefunctions can be defined 
as Xn=Qn[P{z)e-T\^\'\l)). 

Using these expressions it is possible to obtain, following the procedure shown in the case of 
Xi, the appropriate many-body wavefunctions for every value of B and q. In particular all the 
antisymmetric states given by 

N 

j 

(with A^^^ an odd Laughlin state ([8])) are fermionic ground states. Besides, for even values m > 
2n, also the symmetric bosonic wavefunction \E'n"^^ describes a possible ground state for both 
intra-species and inter-species interactions. 

Degeneracy points and non-abelian anyons: An important property of the spectrum Q is that 
there are points in which degeneracies of the Landau levels arise: for g^/ B = l+2n the two lowest 
energy levels cross and the ground state degeneracy of the single particle is doubled. Notice that 
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FIG. 2: Energies of Xn+i (black solid line) and Xn (^^^ dashed line) for n = 0, . . . , 10 as a function of 
/B. The crossings of the ground states happen in the degeneracy points /B = (1 + 2n), denoted by 
solid circles. Inset: Degeneracy points of Landau levels. 



there is complete sequence of points in which two Landau levels have the same energy (integer in 
unit of B) as depicted in the inset of Fig. |2l 

For / B = l + 2n, the single-particle ground states are all the superpositions of wavefunctions 
in Xn Xn+i- Let us consider for simplicity the case of = 3B at the crossing between the 
states Xi and X2 (but our conclusions can be easily generalized to all the other degenerate points). 
Once we introduce the intra-species interaction between atoms, this double degeneracy implies a 
novel form of the multiparticle ground state which is quite different from ([91) : in the degeneracy 
points the antisymmetrization hides a clustering of the particles into two sets of particles A and 
B that are physically different and refer to states Xi X2 ■ For fermions, we can rewrite the 
ground state ^Ic for 2N particles making explicit this clustering: 



n^i^ n -%)n^2;' n 



1 Zi Zj 



2N 



III - i) 



(12) 



.k£A i<jeA leB i<jeB 

ric is obtained through a Slater determinant involving all the states Xi and X2 up to the power 
create a quasi-hole excitation, the total angular momentum must be increased by N/2 
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and, given 2M quasi-holes of coordinates (h, the corresponding wavefunction can be expressed as 

M 



n 



qh 



A 



n ^1 n ^2 n - ^ ) n n - 

.ABA A h=l 



2M 



2JV 

e i;...;) (13) 



B B h=M+l 

These quasi-holes obey a fermionic statistics once two of them in the same cluster (A or B) are 
exchanged; however they show the same fusion rules of the Ising model with defined fermionic 



parity 



2M- We notice that if 



[|23l1 characterizing the Moore and Read (MR) Pfaffian state [H 
one is allowed to consider their linear superpositions at the degeneracy points of this ultracold 
atoms setup (corresponding to atoms in the space spanned by Xi and X2 moving in higher angular 
momentum states), then such superpositions obey non-Abelian braiding rules. 

The state is the ground state corresponding to the highest possible density in the degeneracy 
point given by = 3B; decreasing the filling factor (i.e. the density) the ground states are 
characterized by polynomials with higher degrees; thus it is possible to obtain other ground states 
of the system with 2N particles different from ^ and (fT2l) . but nevertheless characterized by a 
clustering of the particles in Xi and X2 ■ We can introduce, for instance, the wavefunctions: 

1 



^Hf = iGl;iG2;j — Q2;iGl-j 



A 



(m) 
2N 



(14) 



where m is an odd integer, A2'^^ is the Laughlin state ([8]) for 2N particles, and Hf indicates the 
Haffnian, which is a symmetric version of the Pfaffian. The state is antisymmetric over all 
the particles and this guarantees that the inter-species interaction gives a zero contribution, even 
if the wavefunction fijjf doesn't vanish for Zi — )■ Zj, that makes it closely related to the MR state. 
Moreover, for every m > 4 (also for even m), is a ground state for every two-body repulsive 
contact interaction, therefore it is stable under (repulsive) inter-species perturbations. 

Another interesting ground state of T-L (with higher angular momentum) is a deformed MR state 
that can be described by the wavefunction: 

TV 2N 



n 



MR 



5 n^i* n ^2;^ 



pf 



Zi Zj 



A 



(2) 
2N 



(15) 



^i=l i=l+N 

where S is the symmetrization over all the particles and Pf is the Pfaffian operator. This state 



shares all the main characteristics of the MR wavefunction ll25ll . and, in particular, its excitations, 



introduced through a clustering similar to equation (fT3l) . are non-abelian Ising anyons. Notice that 
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the state (fTSl) is formally identical to a state considered in n25\\ . where it was shown that it has 



the same statistical properties of the MR state. In fact, the relation between ^mr and the usual 



MR state 124, 



251 ] allows to relate their norms; therefore, the exchange statistics characterizing 



excitations ( of the kind 



^MR (Ca, Cb, CcXd) = S [\ Qi-i Q2;i 



( - Ca) {Zj - Cb) {Zj - Cc) {Zj -Cd)+i^J \ ^(2) 



Zi Zj 



2N 



is described by four Ising anyons as in the case of the MR state, because the monodromy and Berry 
phases acquired in the exchange of two excitations coincide. 

Conclusions: In this paper we studied ultracold atoms in an artificial non-abelian gauge poten- 
tial: we computed the single-particle energy levels, showing that the non-abelian part of the vector 
potential splits the degeneracy of Landau levels. In presence of strong intra-species interactions, 
deformed Laughlin states and abelian excitations are obtained. However, at the points in which 
the Landau levels are again degenerate, multiparticle ground states characterized by a clustering 
similar to the one identified in Moore and Read states emerge. The system is then suitable to 
have non-abelian anyonic excitations, as the Ising anyons, opening the possibility to tune artificial 
non-abelian gauge potentials in ultracold systems to induce and manipulate non-abelian anyons. 

Given the importance of non-abelian anyons for topological quantum computation, an impor- 
tant question for future work is the study of the braiding rules for the non-abelian anyons obtained 
at the degeneracy points: from this point of view, one could consider more components (e.g., with 
tetrapod schemes) in order to explore the possibility to have lines or eventually regions of de- 



generacy. Similarly interesting would be the characterization of the topological order [|26ll in the 
transitions between deformed Laughlin states and deformed Moore-Read states at the degeneracy 
points, which could show similarities with transitions between Halperin and Moore-Read states in 
bilayer systems 12711 . 
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Appendix A: Artificial non-abelian gauge potential in a rotating tripod system 

The aim of this Appendix is to show how it is possible to obtain the following 2x2 non-abelian 
potential 

B B 
Ax = qcTx - —yl , Ay = qay + —xl (Al) 

starting from a fine tuned system of tripod atoms of the kind described in [fiol (here cr^ and cFy 
are the usual Pauli matrices and I denotes the identity matrix). In particular we will obtain the 
pote ntial A(q) starting from a particular family of configurations of the Rabi frequencies fi^ in 



ttloU 



ni = n sin (9) COS ((/)) e^^i , n2 = n sin (0) sin e'^^ , fig = i7 cos (9) e'^' (A2) 

where 5*1 and 5*2 are functions of the position and of the parameter q while the angles (p and 9 and 
5*3 are chosen constants. 

The key elements to obtain the desired vector potential A (lAll) are the rotation of the whole 



system (in a way similar to the rotating system in [l lql ) and a suitable gauge transformation. Let us 



consider a system of tripod atoms in an inertial frame of reference characterized by a non-abelian 
gauge potential A, a scalar potential Vj-ot = ^{A) + V as the one described in [10] and a harmonic 
confining potential u;r^/4 (here and in the following we will consider m = 1/2 and ^ = 1). The 
corresponding Hamiltonian reads: 

HiF = + i) ' + ^u; + V,ot (A3) 

If we put the whole system in rotation with angular velocity fi the Hamiltonian in the rotating 
frame of reference reads: 

^Rot = (p + i) ' + ^wV^ + Qh, + Kot (A4) 
where we introduced the gauge invariant angular momentum 

L = f X (p + a) (A5) 



and all the coordinates are now considered in the rotating frame. We can now rewrite i^Rot intro- 
ducing the gauge potential: 

~ B B 
Ax = Ax —y , Ay = Ay + —X (A6) 
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where we put B = co. We obtain: 

Hnot = {p + Af + AL, + Kot (A7) 

with A = f2 — cj. 

Our aim is to identify the correct family of Rabi frequencies, V^ot and gauge transformation 
such that: 

H={p + Af + AL, (A8) 

with A given by (|AH) and = rx p being the usual angular momentum in the rotating frame. In 
particular we need: 

A = iqcx^,qay) (A9) 
Vrot = qA{y(j^ - x(jy) (A 10) 

In order to obtain, in the rotating frame, the potential A in (|A1|) starting from the potential A 
given in [|iqI] 

All = cos^ + sin^ 0V5i3 (All) 

A22 = cos^ ^ (cos^ 0V^i3 + sin^ 0V^23) (A12) 
A12 = cose Qsin20V5i2 - zV0^ (A13) 

we need a suitable unitary gauge transformation 0{r). In particular the field transforms as .4 
OAO"^ — iOVO^ and thus we must have 

OAO^ - iOVO^ = A= {qa^, qay) (A14) 

From the definition of A it is evident that, choosing a constant (p, it is not possible to obtain 
Ay oc ay but it is easy to check that we can obtain Ay = kl — qa^ and Ax = qcr^ for a suitable 
choice of the parameters as functions of q. Therefore the gauge transformation we will apply is: 

^ ^ , with O = e'^y-'^""^ (A15) 

with k to be defined in the following. In this way we obtain: 

Ax = qcTx ^ Ax = qcJx , Ay = kl- qa^ ^ Ay = qcTy (A16) 
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We must also consider that the scalar potential in /lO^ are affected by O as O (V^ + $) 0\ thus, in 
order to obtain (IA8I) out of (IA7I) we must have: 



0{V + ^)0^ = V,ot = -Ar X A (A17) 

and then: 

V + <^ = qA{y(T^ + xa,) (A18) 

Now we can find the suitable parameters to satisfy (IA16I) and (IA18I) . First of all we impose 
= 7r/4 and = cost. Then, from the definition of A we obtain that: 

{S^ + S2) = (A19) 

cosOd^iSi- S2) = 2q (A20) 

2A; - 2g = dy (Si + S2) (A21) 

2g + 2A; = cos^ Ody {Si + ^2) (A22) 

A possible solution is given by: 

5i = X{x + y) (A23) 

^2 = X{-x + y) (A24) 

with X = q/ cos 6 in order to satisfy the first two equations; we obtain from the last two equations: 

cos^^ - 2cose - 1 = => cos0 = l-y2 (A25) 

Therefore: 

q q 1 + cos^g 2-^/2 

^ = 7; = 1= y k = — — q = ^ q (A26) 

cos^ l-^/2 2cos^ 1-^2 

In this way we found the values of Si, S2,(p,9 in order to obtain the right A after the gauge 

transformation. Let us consider now the scalar potentials; imposing = 7r/4 and cos 6 = 1 — \f2 

we find from [10] and from (IA18I) : 

1/11 + $11 = ^ ^ ^ + — sin^ 9 = qAx (A27) 
V22 + $22 = " : cos^ e + Vs sin^ ^ + cos^ 6 sin^ 6 = -qAx (A28) 





2 


Vi 


+ V2 




2 


Vi 


-V2 


2 



Vi2 + $i2 = \ COS 9 = qAy (A29) 
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The solution is given by: 

A2 



Vi = qAx + XAy - '— sin^ Q (A30) 
V2 = qAx - XAy - — sin^ 9 (A31) 



4 

= kAx - ^q^ (A32) 



with A and k given by (|A26I) . 

In conclusion, we showed that the values of the angles 



— u = arccos 
4 



(1 - (A33) 



and the linear functions S*!, 5*2, Vi, V2, V'i described above, allow to define the right Rabi frequen- 
cies in order to obtain, through the gauge transformation (IA15I) . the Hamiltonian (|A8I) for a rotating 
system. In particular we are able to describe a family of rotating physical systems, characterized 
by q, with the potential A{q) (|AH) that has been analyzed in our paper. 
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